/* Check JER_Matrix::Inverse(). 
   reed_sol... already checked. */
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include "reed_sol-2.h"
#include "jerasure-2.h"
#include "galois-2.h"

using namespace std;

int main( int argc, char** argv )
{
	JER_Matrix *jm;
	JER_Matrix *jm2;
	int i, j;

	int m = atoi( argv[1] );
	int k = atoi( argv[2] );
	int w = atoi( argv[3] );
	jm = reed_sol_extended_vandermonde_matrix( k, k, w );
	jm->Print();
	fprintf( stdout, "----------\n" );
	jm->Inverse();
	jm->Print();

	return 0;
}
